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We use a combination of a many-body model analysis with an " ab initio" band structure calcula- 
tion to derive the temperature dependent electronic quasiparticle structure of the rare-earth metal 
Gadolinium. As a local-moment system Gd is properly represented by the ferromagnetic (multi- 
band) Kondo-lattice model (s-f (d-f) model). The single-particle part of the model-Hamiltonian is 
taken from an augmented spherical wave (ASW) band calculation. The proposed method avoids the 
double counting of relevant interactions by exploiting an exact limiting case of the model and takes 
into account the correct symmetry of atomic orbitals. The "a priori" only weakly correlated 5d 
conduction bands get via interband exchange coupling to the localized 4f levels a distinct tempera- 
ture dependence which explains by a Rudermann-Kittel-Kasuya-Yosida (RKKY) -type mechanism 
the ferromagnetism of Gd. We get a self-consistently derived Curie temperature of 294.1 K and 
a T = 0- moment of 7.71 /Ub, surprisingly close to the experimental values. The striking induced 
temperature-dependence of the 5d conduction bands explains respective photoemission data. The 
only parameter of the theory (interband exchange coupling J) is uniquely fixed by the band calcu- 
lation. 

PACS numbers: 71.10.-w,75.30.Et,71.20.Eh 
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I. INTRODUCTION 



The rare-earth metal Gd is one of the four elemen- 
tal ferromagnetic metals; the others are Fe, Co, Ni. It 
crystallizes in the hep structure with a lattice constant 
a = 3.629 A, c/a — 1.597 1 . Magnetic properties re- 
sult from the half-filled 4f shell (L = 0,J = 5 = |) 
which gives rise to strictly localized magnetic moments. 
Conductivity properties are due to partially filled 5d/6s 
conduction bands. As to the purely magnetic properties 
Gd is considered an almost ideal Heisenberg ferromag- 
net with a Curie temperature of Tq = 293.2 K and a 
zero-temperature moment of /j,(T = 0) = 7.63 /xb^- The 
latter indicates an induced polarization of the conduction 
bands of at least 0.63 (j,b due to an interband exchange 
coupling between itinerant band electrons and localized 
4f electrons. The rather strict localization of the 4f wave 
function^ prevents a sufficient overlap for a direct ex- 
change interaction between the 4f moments. The cou- 
pling between the moments is therefore mediated by po- 
larized 5d/6s conduction electrons (RKKY), i.e. strongly 
influenced by the electronic structure. 

Although the ferromagnetic ground-state of Gd is of 
course without any doubt, it is still a matter of debate 
how to get this fact by an 11 ab initio" band structure 
calculation. Numerous investigations of the electronic 
ground-state properties of Gd have been performed in the 
recent p^^^^^iiSilLH, a ll i n the framework of den- 
sity functional theory (DFT). They provide a convincing 
description of ground-state properties such as the lattice 
constant, the hep-crystal stability, the c/a ratio, the mag- 
netic moment^, the bulk modulus, and the Fermi sur- 
face paramtersiS. On the other hand, a standard local 



density approach (LDA) to DFT predicts an antiferro- 
magnetic ground state if the 4f electrons are considered 
as valence electrons. In a detailed analysis Kurz et aliii 
have demonstrated that the reason for the incorrect pre- 
diction of antiferromagnetism is the wellknown difficulty 
of LDA correctly to describe strongly localized electrons. 
The LDA calculation of Ref. ljj poses the nearly disper- 
sionsless majority 4f bands some 4.5 eV below the Fermi 
energy while the minority 4f bands are directly above the 
Fermi energy leading to a certain itinerancy of the 4f elec- 
trons. These findings are at variance with the results of 
combined direct (XPS) and inverse (BIS) photoemission 
experiments^ which observe occupied 4f f-states at the 
binding energy (—7.44 ± 0.1) eV and unoccupied 4f in- 
states at (+4.04 ± 0.2) eV, i. e. distinctly away from the 
Fermi edge. A special consequence of the wrong position 
of the down spin 4f states is an extremely high density of 
states close to the Fermi edge and therewith an unrealis- 
tically big 7-value of the electronic heat capacity. The 
most important consequence of the wrong 4f |-position, 
however, is the prediction of antiferromagnetism. All cal- 
culations, which treat the 4f electrons as valence elec- 
trons, irrespective of whether the LDA or the general- 
ized gradient approximation (GGA) is applied, end up 
with an antiferromagnetic Gd-ground-stat o 11 : 14 . Today 
it is clear how to remove this inadequacy of LDA (GGA) 
"by hand"—. One has to remove the 4f |-states from 
the Fermi energy. This can be done simply by consider- 
ing the 4f states as "core states", so that they are not 
allowed to hybridize with any other states on neighbour- 
ing atomsii* 1 ^. Another way is to apply the so-called 
LDA+C/ methodic, which introduces strong intraatomic 
interactions of the localized states in a Hartree-Fock-like 
manner. The main effect is a splitting apart of the occu- 
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pied and unoccupied 4f statesi2iiLi&, i.e. in particular a 
removal of the minority 4f states from the Fermi edge. 

Needles to say that in a real "ab initio" DFT-procedure 
all electron states, i.e. in particular the 4f states, should 
be treated as valence states. To declare the 4f states 
as core-states or to introduce at a convenient place the 
"Hubbard-U " surely corrumpes a bit the "first princi- 
ples'''' character of the band calculation. The only mo- 
tivation is to compensate the LDA artifact which pre- 
vents the correct ferromagnetic Gd-ground state. We 
have recalculated the Gd-band structure using a new im- 
plementation of the augmented spherical wave (ASW) 
methodize. The 4f electrons have been treated as va- 
lence electrons and the scalar-relativistic approximation 
of the Dirac equation has been used. Although the an- 
tiferromagnetic configuration turns out to minimize the 
total energy, the ferromagnetic order has been assumed, 
firstly because it is closer to reality, secondly because we 
need these data for our further procedure. The results 
for some high-symmetry directions are plotted in Fig. ^ 
We recognize the wrong position of the rather flat 4f dis- 
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FIG. 1: Spin resolved T = band structure of ferromagnetic 
Gd as a function of the wave vector, obtained by a scalar- 
relativistic ASW calculation. Solid lines for ] states, dotted 
lines for [ states. The energy zero is defined by the Fermi 
edge. The flat dispersions are the 4f levels. 

persions. The fairly broad 5d/6s dispersions exhibit an 
exchange splitting with a weak k-dependence. The re- 
sult is in good agreement with other first-principles cal- 
culations that treat the 4f states as valence states&Siii. 
The occupied and almost dispersionless 4f f-bands pro- 
vide the major part of the magnetic moment (w 7 /xb) 
while the 5d splitting at the Fermi edge accounts for the 
experimentally observed excess moment of w 0.63 ^b- 

The mentioned 5d/6s exchange splitting must be in- 
duced since the 5d/6s electrons can be considered as 
only weakly correlated and "a priori" non-magnetic. The 
splitting is obviously due to a strong interband exchange 
interaction between the flat 4f states and the extended 
5d/6s conduction states. Precondition for that is a ferro- 
magnetic order of the localized magnetic moments built 



by the half-filled 4f shells. The next neighbour distance is 
too large for a direct exchange interaction. The moment 
coupling is of indirect nature mediated by a spinpolar- 
ization of the conduction electrons due to the mentioned 
interband exchange interaction with the localized 4f elec- 
trons. This so-called RKKY mechanism is strongly de- 
pending on the electronic structure. To understand the 
ferromagnetism of Gd does therefore mean first of all to 
understand its electronic structure. 

The induced 5d/6s exchange splitting is still a mat- 
ter of controversial debate, in particular what concerns 
its temperature dependence^. Is it collapsing or non- 
collapsing for T — > To? Photoemission data appear to 
be not unique. Some experiments point to a collaps- 
ing ("Stoner-like") behavior—, others exhibit a splitting 
that does not shift very much with temperature ("spin- 
mixing") persisting in the paramagnetic phase 20 . In the 
latter case the demagnetization for T —* Tq is reached 
by a redistribution of spectral weight rather than by a 
gradually increasing overlap of respective spin peaks. To 
find out what is really going on one needs a theory for the 
full temperature-dependence of the electronic structure. 
Pure " ab initio" band calculations are restricted to T = 
being therefore insufficient for this purpose. Sandratskii 
and Kubler— have proposed a DFT-based theory where 
finite-temperature effects are simulated to a certain de- 
gree by a respective directional disorder of the spatially 
localized 4f moments. Even though being an interest- 
ing ansatz it certainly cannot replace the full statistical 
mechanics of the local-moment ferromagnet. 

A key-quantity of ferromagnetism is the Curie tem- 
perature Tq, It is the aim of each theory for a ferromag- 
netic material to approach Tq as quantitatively as pos- 
sible. On the other hand, it is a very sensitive term to 
get. Several attempts have been started to estimate Tq 
from total-energy calculations by use of the LDA-DFT 
schemeii* 2 ^^. For this purpose the energy data are in- 
serted into simple mean-field formula for the magnetic 
transition temperature, very often arriving at astonish- 
ingly accurate Tc-values. However, it is surely not unfair 
to state that such estimates cannot replace a full theory 
of the Gd ferromagnetism. The latter requires access to 
an electronic structure calculation which fully accounts 
for decisive temperature and correlations effects. To our 
information such a complete theory does not yet exist for 
the prototype local-moment ferromagnet Gd. It is the 
aim of this paper to present a methodical approach to 
the temperature-dependent electronic structure of ferro- 
magnetic local-moment metals with a direct application 
to Gd. 

We present a theory of the electronic quasiparticle 
structure of the ferromagnetic 4f metal Gd that yields 
in a self-consistent manner the electronic as well as the 
magnetic properties. The approach shall work for arbi- 
trary temperatures regarding in particular electron cor- 
relations effects. To get in this sense a realistic picture of 
Gd we combine a "first principles" band structure calcu- 
lation with a many-body evaluation of a properly chosen 
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theoretical model similar to previous work on the band 
ferromagnets Fe, Co, and Ni24. We consider the underly- 
ing proposal a continuation and extension of a previously 
pubished paper™ that already dealt with the electronic 
quasiparticle structure of Gd. However, in the previous 
case we did not succeed in getting the ferromagnctism 
self-consistently, i.e. simply via the special electronic 
structure. There appeared a serious ambiguity how to 
handle the d-band degeneracy, i.e. how to perform the 
necessary decomposition of 5d-band in non-degenerate 
subbands. The decomposition in Ref. |25| did not em- 
phasize the correct symmetry of atomic orbitals. Sepa- 
rate model calculation a 26 ! 27 revealed that the then-used 
d-band decomposition is inconvenient for a ferromagnetic 
order of the local 4f moments. We propose in this paper 
a new ansatz by which one gets correctly the electronic 
structure as well as the magnetic order of Gd. 

The general procedure is briefly described in the next 
section. Central part of the procedure is a many-body 
evaluation of a properly chosen theoretical model. In sec- 
tion IIIII we introduce and justify the (multiband) Kondo- 
lattice model (KLM) as a good starting point for an 
at least qualitative understanding of local-moment ferro- 
magnets such as Gd. We explain how to combine it with 
an LDA-DFT band calculation to come to quantitative 
statements. The KLM provokes a non-trivial many-body 
problem which for the general case cannot be treated 
rigorously. In section llVl our theoretical approach is rep- 
resented. In the last step (section 0) we combine the 
model analysis with a band structure calculation to get 
the electronic quasiparticle spectrum of Gd and its tem- 
perature dependence, that, on the other hand, fixes the 
magnetic properties of the rare earth metal as, e.g., the 
Curie temperature and the magnetic moment. 

II. GENERAL PROCEDURE 

Our study aims at a quantitative determination of 
the temperature-dependent electronic structure of the 4f- 
ferromagnet Gd. The general concept is rather straight- 
forward and consists of three steps. The important first 
step is the choice of a suitable theoretical model. The 
main physics is due to the existence and the mutual influ- 
ence of two well defined subsystems, quasi-free electrons 
in rather broad conductions bands (5d/6s) and localized 
electrons with extremely flat dispersions (4f). The theo- 
retical model is defined by its Hamiltonian; 

H = H + H 1 (1) 

more strictly, by its interaction part Hi . This particular 
operator shall incorporate all those interactions which 
are responsible for the characteristic phenomena under 
study. In the present case Hi should cover the decisive 
electron correlations which determine the magnetic prop- 
erties and the characteristic temperature dependence of 
the electron quasiparticle spectrum. Our proposal for Hi 
is discussed in the next section. 



While there is no contribution of the f electrons to the 
kinetic energy, the part of the band electrons reads: 

ijamra' 

= E ( T k mm ' - ^w)4 mCT c km/(T (2) 

kcrmm' 

c jma( c jma) ^ s the creation (annihilation) operator for a 
Wannier electron at site Rj in the orbital m with spin a 
(a =t, I). c kmcr (c kmcr ) is the respective Fourier transform 

c kmi r— ^ ] c jma e 3 (3) 



The hopping process from site Ry to site Ri may be 
accompanied by an orbital change (m' — > m). T™ m are 
the respective hopping integrals: 

k 

fi is the chemical potential. The single-particle part Ho of 
the model-Hamiltonian H stands, as usually, for the ki- 
netic energy of the itinerant charge carriers and for their 
interaction with the lattice potential. However, it shall 
furthermore account for all those interactions which are 
not directly covered by Hi . By construction these inter- 
actions are not important for the magnetic properties and 
the temperature-dependence of the electronic spectrum 
of the underlying ferromagnetic material (Gd). Never- 
theless they may influence the rough structure of the 
spectrum so that they cannot be neglected if our study 
really aims at a quantitative description of Gd. For this 
reason we perform in the second step of our procedure a 
full selfconsistent band structure calculation within the 
LDA-DFT scheme in order to replace the single-particle 
energies T™ in Hq by the effective energies of the LDA: 

jW > rpmm' (LDA) (6) 

Since the "other" interactions are, by construction of 
the model, not responsible for the temperature effects, 
we can incorporate them for any temperature, e.g. for 
T = where LDA works. It is therefore guaranteed 
that all the other interactions are well accounted by the 
single particle energy renormalization JHJ. However, the 
problem of a double counting of just the relevant interac- 
tions, namely once explicitely in Hi and then once more 
implicitcly in Hq, must carefully be avoided. How we 
circumvent this problem in the case of the local-moment 
ferromagnet Gd is explained at a later stage. 
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In the third step we apply a many-body formalism in 
order to investigate how the effective single-particle en- 
ergies change under the influence of the "relevant" inter- 
action Hi into temperature -, carrier concentration (n) - 
and possibly spin - dependent quasiparticle energies: 

Hi : T? m ' (LD A) — » E ma (k, T, n) (7) 

The quasiparticle energies are derived as the poles of the 
single-electron Green-function matrix: 



G k<7 (E) = h (£ + M )I-T k -M k(7 (£) 



(8) 



M~ka{E) is the selfenergy matrix, the determination of 
which solves the problem. T k is the hopping- matrix. The 
elements of the Green-function matrix are the retarded 
single-electron Green functions: 

r+oa 

(K ma ,c+ m , a )) E = -i / dte* Et ({c kma (t),c+ n , a (0)} + ) 
Jo 

(9) 

[...,.. •]+(-) means the anticommutator (commutator) 
and (...) is the thermodynamic average. The elements 
of the selfenergy matrix formally solve the Green function 
equation of motion: 

(([c W) ffi]_;c+ mV )) E ^Y, M ^ m ''(E)((c km ^;c+ m>a )) E 

m" 

(10) 

We will discuss our results in terms of spectral densi- 
ties (SD) and quasiparticle densities of states (Q-DOS), 
because both have a direct relationship to the experi- 
ment. Except for respective transition matrix elements 
the spectral density expresses the bare line shape of an 
angle- and spin - resolved (direct or inverse) photoemis- 
sion spectrum: 

S km *(E - /j)=--ImG k T(£ - ,i) 

7T 

= -^ Im ((Ckm ff ;cL»B-p (11) 

An additional k-summation yields the quasiparticle den- 
sity of states (angle-averaged photoemission spectrum !) 



(12) 



that in general will be temperature -, carrier concentra- 
tion -, lattice structure -, and in particular for ferromag- 
netic systems, explicitely spin-dependent. 



III. THEORETICAL MODEL 

A. Model-Hamiltonian 

We still have to fix the interaction part Hi of the model 
Hamiltonian ||TJ for the correlated system of localized (4f ) 



and delocalized (5d) electrons. We presume from the 
very beginning an on-site Coulomb interaction between 
electrons of different subbands, 



L\ ■ ■ -L/4 <7(7 



U Li-L 4 cl ia c+ 2a ,c L3CT ,c Lia (13) 



For simplicity we drop for the moment the site index i 
and the respective summation, that will be reintroduced 
at the end of the following consideration. L\ . . .L\ de- 
note the different bands, and Ul 1 ---l 4 , are the Coulomb 
matrix elements. Restricting the electron scattering pro- 
cesses caused by the Coulomb interaction to two involved 
subbands, only, we get instead of (|13fl : 



Hi 



1 



'■La C L'(7' C L'a' C La 



LL' aa f 

+ JLL'C La C L , (y ,C L(J ,C L , (J 

In the case of Gd the band indices L and L' can be at- 
tributed either to a flat 4f band (L — > /) or to a broad 
(5d/6s) conduction band (L — > to). In an obvious man- 
ner we can then split the Coulomb interaction into three 
different parts, 



Hi = H, 



dd ' 



H 



ff 



H 



(15) 



depending on wether both interacting particles stem from 
a conduction band, Hdd, or both from a flat band -ff//, or 
one from a flat band the other from a conduction band, 
Hdf- The first term, Hdd, refers to electron correlations in 
the broad conduction bands. We consider them not to be 
decisive for the characteristic Gd-physics. According to 
our concept (section[nJ Hdd does not enter explicitely our 
model being rather accounted for by the single-particle 
energy renormalization ©. Hff is built by pure 4f cor- 
relations. The main influence of the 4f electrons on the 
Gd-physics is due to the fact that they form permanent 
localized magnetic moments. So Hff is unimportant as 
part of our model-Hamiltonian and we are left with the 
interaction between localized and itinerant electrons: 

Hdf = y ] \P mfC^naC f G ' c f a' C ma + Jmf 'Cma c f a' C m<r' c / a 
mfaa' 

+ 2^rnf c m<j C mcr' C fcr' C fcr + 2 ^fm C ~fcr C ~fa' C m<j' c ma } (16) 

The last two terms do not contribute since the Gd 3+ -4f 
shell has its maximum spin S = 7/2. All the seven 4f 
electrons have to occupy different subbands and none of 
the seven subbands will be doubly occupied. By use of 
the electron spin operator, 

cr+ = hc^c^ cr" = ?ic^c T ; a z = |(n T - nj (17) 

(n a = c^c a ) we get Hdf in the following compact form 

H d f = -^2^ J m f(Tm ' & f + (U m f - -J m f)n m n f 



mf 



»i f 



(18) 
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with n TO (/) = n. m (/)f + n m (m. For all processes of in- 
terest the number of f electrons per site is fixed, rif is 
therefore only a c-number. The last term in 1)18(1 does not 
really provide an fd-interaction. It leads only to a rigid 
shift of the atomic levels being therefore fully accounted 
for by the renormalization © of the single-particle part 
of the Hamiltonian. By defining the spin operator S of 
the local f moment 



7> 



(19) 



and by assuming that the interband exchange J m f is in- 
dependent of the special index-pair to, / 



(20) 



the interaction term reads after the reintroduction of the 
lattice site dependence: 

J 



Hdf — ~2 ^ ' & irn ' 



J 

2h 



iin — atma 



Here we have used the abbreviations 



5? 



} (21) 



(22) 



The single-band version (non-degenerate s band) of 1(21(1 
is well-known as the interaction part of the so-called 
Kondo-lattice model (KLM) 28 , in the older literature 
more appropriately denoted as s-f or s-d mode&22£&. In 
the multiband case we have in Hdf simply an additional 
summation over the orbital index to. The first term of 
((21(1 describes an Ising-like interaction of the two spin 
operators, while the other provides spin exchange pro- 
cesses between localized moment and itinerant electron. 
Spin exchange may happen by three different elementary 
processes: Magnon emission by an itinerant .[.-electron, 
magnon absorption by a t -electron and also formation 
of a quasiparticle, which is called "magnetic polarori'. 
The latter can be understood as a propagating electron 
"dressed" by a virtual cloud of repeatedly emitted and 
reabsorbed magnons corresponding to a polarization of 
the immediate localized spin neighbourhood. 

Our model-Hamiltonian, built up by the partial oper- 
ators @ and 



H = H Q + H, 



df 



(23) 



can be considered as "multiband Kondo-lattice model" 
(m-KLM). While in the interaction part Hdf the multi- 
band aspect appears only as an additional summation, 
the subbands are intercorrelated via the single-particle 
term H . 

An important model parameter is of course the effec- 
tive coupling constant JS/W where W is the width of the 



"free" Bloch-band and S the local spin value. It turns 
out that in particular the sign of J is decisive. Other 
model parameters are the lattice structure and the band 
occupation 

n = J2( n rna) (24) 
ma 

In case of an s-band n is a number in between and 2. 

B. Exact limiting case 

The many-body problem provoked by the model 
Hamiltonian ((23(1 is rather sophisticated, up to now not 
exactly solvable for the general case. Fortunately, how- 
ever, there exists a non-trivial, very illustrative limit- 
ing case which is rigorously tractable, nevertheless ex- 
hibiting all the above mentioned elementary excitations 
processes2ii22ii£. It refers to a single electron in an oth- 
erwise empty conduction band being coupled to a ferro- 
magnetically saturated moment system. Such a situation 
is met, e.g., for the ferromagnetic semiconductor EuO at 
T = 0. Because of the empty band and the totally aligned 
spin system the hierarchy of equation of motions of the 
single-electron Green function (JSJ decouples exactly. 
One can exploit exact relationships of the following kind: 

(. • ■ c ma ) = (c+ na ...)= 0; (...#) = ...) = 
(...£?>=<£?...)= ftS (...) (25) 

A troublesome but straightforward calculation then ar- 
rives at the following result for the selfenergy matrix © 



M ka (E) = --z a JSi+(l-z a ) 



\J 2 S j-G a (E+±JS) 

I-\J±G {E+\JS) 
(26) 



-Go(E) 



(27) 



The | spectrum is especially simple because the f elec- 
tron cannot exchange its spin with the parallely aligned 
local spin system. Only the Ising-type interaction in 1(21(1 
takes care for a rigid shift of the selfenergy by —^JS. The 
spectral densities l(TT|) arc ^-functions representing quasi- 
particles with infinite lifetimes. Real correlation effects 
appear, however, in the J, spectrum. The essentials can 
be seen already for a non-degenerate s band. Fig.[21shows 
the energy dependence of the |-spectral density Ski C^O 
for some symmetry points. Furthermore, we have chosen 
a sc lattice, S = 1/2 and W = 1 eV. 

For weak coupling (e.g. J = 0.05 eV) the spectral 
density consists of a single pronounced peak. The fi- 
nite width points to a finite quasiparticle lifetime due to 
some spin flip processes, but the sharpness of the peaks 
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indicates a long living quasiparticle. This changes dras- 
tically even for rather moderate effective exchange cou- 
plings JS/W. One observes in certain parts of the Bril- 
louin zone, for strongly coupled systems even in the whole 
Brillouin zone, that the spectral density splits into two 
parts. The sharp high-energy peak belongs to the forma- 
tion of the magnetic polaron while the broad low-energy 
part consists of scattering states due to magnon emis- 
sion by the J. electron. As long as the polaron peak is 
above the scattering spectrum the quasiparticle has even 
an infinite lifetime. The scattering spectrum is in gen- 
eral rather broad because the emitted magnon can carry 
away any wave-vector from the first Brillouin zone. Be- 
cause of the concomitant spinflip magnon emission can 
happen only if there are | states within reach. Therefore, 
the scattering part extends just over that energy region 
where P](E) ^ 0. Sometimes, as e.g. for J = 0.6 eV 
at the r point (Fig. the scattering part is surpris- 
ingly bunched together to a prominent peak, therefore 
certainly visible in a respective (inverse) photoemission 
experiment. Note that the results in Fig. |3 are exact and 
free of any uncontrollable approximation. They exhibit 
typical correlation effects which are by no means repro- 
duceable by a single-electron theory. 

Very important for the following procedure is the sim- 
ple | result. It tells us that at (T = 0, n = 0) the 
df-exchange interaction takes care only for a rigid shift 
of the total energy spectrum without any deformation, 
which is therefore identical to the "free" Bloch spectrum. 
Furthermore for this special case a mean-field approxima- 
tion turns out to be exact. Though not exactly provable, 
many reliable approaches^^*^ show that this holds, at 
least to a good approximation, for finite band occupa- 
tion, too. This will be demonstrated in Fig. |3| for the 
actual case of Gd. 



C. Band-structure calculations 

As described in section ITfl the hopping integrals T^ m 
in the single-particle Hamiltonian Ho (J2J) have to ab- 
sorb the influences of all those interactions which are 
not directly covered by our model Hamiltonian (|23|) . For 
this purpose we have performed a spin polarized scalar- 
relativistic ASW-band calculation for ferromagnetic hep 
Gd. The result is plotted in Fig. ^ The 4f electrons are 
considered as valence electrons. In section [I] we already 
commented on the wrong position of the 4f levels. How- 
ever, we need as input for the determination of the quasi- 
particle spectrum only the | part of the bandstructure be- 
cause of the following reason. The main problem, when 
using LDA-DFT results as renormalized single-particle 
input, is to avoid a double counting of just the decisive 
interband exchange coupling, namely once explicitely in 
the interaction-part Hi of the modcl-Hamiltonian and 
then once more implicitely by the effective single-particle 
input. The exact limiting case of the last section gives 
the hint how to circumvent this double-counting problem. 
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FIG. 2: Exact ! spectral density of the single-band Kondo- 
lattice model at T = as function of the energy for several 
symmetry points in the first Brillouin zone and for different 
exchange couplings J. Parameters: S = 1/2, W — 1 eV, 
n = 0, sc lattice. 



For an empty band (n = 0) and ferromagnetic satura- 
tion the | spectrum is only rigidly shifted by a constant 
energy amount (—^JS) compared to the "free" Bloch 
spectrum. As mentioned, model approaches convincingly 
demonstrate that this remains true for finite band occu- 
pation (less than half-filled bands !), too, at least to a 
very good approximation. For less than half-filled bands 
we therefore can identify the (T = 0, a =f)-LDA results 
(without the 4f part) with the free Bloch energies. By 
this procedure we do not explicitely switch off the df- 
interband exchange but rather exploit the fact that for 
the mentioned special case the interband exchange leads 
only to a trivial rigid shift. The wrong position of the 
LDA-4f I states therefore does not bother us because we 
need excludingly the f spectrum. On the other hand, we 
can be sure that all the other interactions, as e.g. the 
Coulomb interaction of the 5d electrons, are excellently 
accounted for by the LDA. We note in passing, that for 
more than half-filled bands we would have to take the 
(T = 0, a =|)-spectrum (particle- hole symmetry). 

After defining the single-particle input there remains 
only one parameter, namely the exchange coupling J. It 
is not considered as a free parameter, but taken from the 
bandstructure calculation. It is commonly accepted that 
an LDA treatment of ferromagntism is quite compatible 
with a mean-field ansataiZ^, so that the exchange split- 
ting in Fig. ^ should amount to A = JS (see next sec- 
tion). We realize, however, that the assumption of a rigid 
splitting is too simple. A slight energy and wave-vector 
dependence of the exchange splitting is found by LDA, 
too. We have therefore averaged the T = 0-splitting over 
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N p prominent features in the Q-DOS of Gd arriving at and the "spin- flip function": 



J = 



N P S 



0.3 eV 



(28) 



This is of the same order of magnitude as found for Eu 
chalcogenides 3 ^!. There are then no other parameters 
in our theory. 

The validity of the above assumptions will be demon- 
strated later with Fig. |3| 



IV. MANY-BODY EVALUATION 

The multiband-KLM l|23[> does not allow a rigorous 
solution except for some special cases. Approximations 
are unavoidable. There are two partial problems to be 
solved, one concernig the ferromagnetism mainly pro- 
voked by the localized 4f moments, the other dealing with 
the temperature reaction of the conduction-band states 
due to the magnetic state of the moment system. In addi- 
tion, both parts have to be combined in a self-consistent 
manner. 

For the electronic part we have developed in the past 
several approache o 26 i 34 i 38 , which all lead, at least qual- 
itatively, to the same result. The "interpolating self en- 
ergy approach" 34 ' 38 is in particular trustworthy for al- 
most empty or almost full bands. For intermediate fill- 
ings as in the case of Gd the "moment conserving de- 
coupling approach" (MCDA) 2 - seems to be more recom- 
mendable. So we use it here. Since this approach has 
been exhibited in detail in Ref. ^ we can restrict our- 
selves in the following to the central parts which are vi- 
tal for the understanding of the underlying procedure. 
To make the representation as clear as possible we de- 
velop the method in the next section for the special case 
of a non-degenerate band. The orbital index is then su- 
perfluous. The generalization for the degenerate case is 
straightforward. The investigation of the magnetic part 
follows in the second subsection. 



A. Conduction electron selfenergy 

According to JSJ the solution of the problem needs 
the knowledge of the selfenergy M-^ a (E). The 
above-mentioned MCDA is a non-perturbational Green- 
function theory. It starts from the equation of motion of 
the site-dependent single-electron Green function (JSJ): 

/X(E + fi)8im - T im )G mjcr {E) = hdij 

m 

--j(z a T iitj(J (E) + F ii>jCT (£)) (29) 

Our approximation attacks the equations of motion of 
the "Ising-function" 



F im j a (E) - ((Si cr c m _ cr ;ct))j 



(31) 



These equations of motion contain still higher Green 
functions which are decoupled to get a closed system of 
equations. Let us exemplify the procedure by a "higher" 
Green function of the type ((AJc^, i?<^]_; c^)) e, where 
A i is any combination of local-moment and band oper- 
ators. The off-diagonal terms i ^ m are approached by 
use of the selfenergy elements Mi ra ( "selfenergy trick"), 
in a certain sense as a generalization of the exact equation 



((A i [c la ,H df Uc+)) 



M lra (E)((A lCra -c\ a )) 



(32) 

The right-hand side is a linear combination of "lower" 
Green functions with the selfenergy elements as selfcon- 
sistently to be determined coefficients. To account for 
the strong local correlations the diagonal terms i — I are 
handled with special care: 

((Ai [c la ,H df ]^; ct a )) E =a (T G ija (E) + (3 a T Hja (E) 

+^F ii>ja (E) (33) 

Such an ansatz is constructed in such a way that all 
known exact limiting cases (atomic limit, ferromagnetic 
saturation, local spin S = 1/2, n = 0, n — 2,. . . ) are ex- 
actly fulfilled. The at first unknown coefficients a a , j3 a , 
7 CT are eventually found by equating exact high-energy 
expansion (spectral moments) of the selfenergy. As the 
other above-mentioned methods^ 3 ^ the MCDA arrives 
at the following structure of the selfenergy: 



M k(T (E) 



1 



Jz a (S z ) + J 2 D^(E;J) 



(34) 



Fi m ,ja{E) — {{Si c ma ; Cj a )) 



(30) 



Restriction to the first term, only, is just the mean-field 
approach to the KLM, which is correct for sufficiently 
weak couplings J, being mainly due to the Ising-part in 
(|21|l . Without the second part it would give rise to a spin- 
polarized splitting of the conduction band. The term 
D)za{E] J) is more complicated being predominantly de- 
termined by spin exchange processes due to the spin-flip 
term in the Hamiltonian (|21|l . It is a complicated func- 
tional of the selfenergy itself, and that for both spin di- 
rections, i.e. (|34|l is an implicit equation for M^ (J (E) 
and not at all an analytic solution. Dk a (E; J) depends, 
furthermore, on mixed spin correlations such as (Sfnia), 
(Sf c^c^) ,. . . , built up by combinations of localized-spin 
and itinerant-electron operators. Fortunately, all these 
mixed correlations can rigorously be expressed via the 
spectral theorem by any of the Green functions involved 
in the hierarchy of the MCDA. However, there are also 
pure local-moment correlation functions of the form (Sf), 
(SfS^ 1 ), ((Sf) 2 ),. . .which also have to be expressed by 
the electronic selfenergy M-^ a (E). 
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B. Modified RKKY interaction 

To get such expectation values of local-spin combina- 
tions we map the interband exchange operator (|21|l on 
an effective Heisenberg Hamiltoniar. 2 >Si2i: 



H d f = 



J 



-2 / j 

ira 



■Si 



elf 



S 7 ' So 



(35) 



We use here again the full multiband version. The map- 
ping is done by averaging out the electronic degrees of 
freedom <7i m — > ((Timr' ■ That means, in the last analy- 



"k+qmcr^kmtr'/ 



\(c) 



sis, to determine the expectation value 
The averaging (■ • • )( c ) has to be done in the conduction 
electron subspace where the local spins Si can be treated 
as classical variables: 

1 



( C k+qmtT C kn 



F7 Tr ( e 



-0H' + 



'k+qmu kmir 



) (36) 



H' is formally the same as in Iffiljl. except for the fact 
that for the averaging process the f-spin operators are to 
be considered as c numbers, therefore not affecting the 
trace. 5' is the corresponding grand partition function. 
We use the spectral theorem for the "restricted" Green 
function, 



Gram 
ka 



(( c kmcr'i C k+qm'cr)) E 



(37) 



,k+q<ri E ) 

to fix the expectation value l|36|l . Eq. I|37|l stands for the 
usual definition |(?JJ of a retarded Green function, only the 
averages have to be done in the Hilbert space of H'. The 
equation of motion of G reads (in matrix representation 
with respect to the orbital indices m, m'): 



G 



(E) 



ka - ' ,k+q(T 

ia"k' 
_ p4 (k'-(k+q))-R I 



(0) 



5q,0<W'G k 



(E) 



i(k-k')-R; 



■Gl 0) (E)G k , 



a' ,k+qcr 



(E) 



(Si ■ a )o-" CT (?kcr',k'. 



(38) 



This equation is exact and can be iterated up to any 
desired accuracy. G k ° (E) is the Green function matrix 
of the "free" electron system: 



G ( °\E) = h (£ + /i)I-T k 



(39) 



If we stop the iteration in (|38l) after the first nontriv- 
ial step, i.e. replacing G on the right-hand side by the 
"free" Green function matrix, then we arrive at the well- 
known RKKY-resuh»2Z., which can be equivalently derived 
by use of conventional second-order perturbation theory 
with respect to J starting from the unpolarized conduc- 
tion electron gas. To incorporate the exchange-induced 
conduction electron spin polarization to a higher degree 
we replace the restricted Green function on the right- 
hand side of not by the "free" but by the full single- 
electron Green function matrix G-k a (E) defined in JHJ: 

Gk'a",k+qa{E) ><5k' ,k+qO~cr" aGk+q<j (E) (40) 

Gka' .k' a" (.E)— W' Gka' (E) (41) 



After some manipulations that replacement leads to the 
following effective exchange integrals: 



T 2 

jctf_ J 

ij 8Ntt 



/+oo 
dEf_(E)* 
-oo 



kqmcr 



Im[(G k .(^-/i)4+ q (^-M) 
G< D) (£-^(£-ri 



(42) 



These effective exchange integrals are functionals of the 
electronic selfenergy M-^ a (E) getting therewith a distinct 
temperature - and carrier concentration dependence. Ne- 
glecting Mka(E), i.e. replacing in Q42J1 the full by the 
"free" Green function, leads to the multiband version 
of the conventional RKKY-exchange integrals 39,40 . Via 
Mka(E) higher order terms of the conduction electron 
spin polarization enter the "modified" RKKY l|42(l which 
is therefore not restricted to weak couplings, only. 

To get from the effective Heisenberg Hamiltonian H35JI 
the magnetic properties of the multiband KLM we apply 
the standard Tyablikow-approximation 41 which is known 
to yield convincing results in the low as well as high tem- 
perature region. All the above mentioned local-moment 
correlations are then expressed by the electronic selfen- 
ergy. We therefore end up with a closed system of equa- 
tions that can be solved self-consistently for all quantities 
of interest. For a detailed discussion of the so- found prop- 
erties of the single-band KLM the reader is referred to 
our previous publications2£*2L2£. We use the theory in 
the next section to find the electronic and the magnetic 
properties of the ferromagnetic 4f-metal Gadolinium. 



V. MAGNETIC AND ELECTRONIC 
PROPERTIES OF GADOLINIUM 

Fig-Elshows the partial (5d, 6s, 6p) quasiparticle den- 
sities of states at T = 0, as they are found by our method 
and compared to the pure ASW-LDA. The | spin parts 
are almost identical for both methods. That confirms 
our procedure, explained in section Till CI for the combi- 
nation of the many-body model evaluation and the "first 
principles" bandstructure calculation. Obviously a dou- 
ble counting of any decisive interaction has almost per- 
fectly been avoided. The still observable very small devi- 
ations might be due to the finite band occupation. The 
statement that the up-spin spectrum at T = is only 
rigidly shifted l|26|l compared to the free spectrum can 
be proven, strictly speaking, only for empty bands. As 
mentioned, a lot of reliable approachesSSiSL 3 ^ support the 
assumption that this is true, at least to a very good ap- 
proximation, for finite carrier densities, too. However, 
slight deviations may appear. Furthermore, the band- 
and wave-vector-independence of the exchange coupling 
J (|28|l is surely an oversimplification and may also con- 
tribute to the deviations in the f spectrum. Nevertheless, 
the almost complete coincidence between LDA and model 
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FIG. 3: Spin resolved densities of states of the 6s, 6p, and 5d 
bands of Gd as functions of the energy at T — 0. Full lines 
for the ASW-LDA calculation, broken lines for our theory, j 
(J.) spectra in upper (lower) halves of the figures. 



results demonstrate that there are hardly any exchange- 
caused correlation effects in the t spectrum of the local- 
moment ferromagnet at T — (ferromagnetic satura- 
tion) . 

The J, part of the T — 0-spectrum, however, exhibits al- 
ready strong correlation effects due to the exchange cou- 
pling of the band states to the 4f moment system, pre- 
dominantly in the 5d subband. They follow from magnon 
emission processes of the down-spin electrons and to a 
lesser extent from the formation of magnetic polarons. 

Integration up to the Fermi edge yields the T = 
contribution of the conduction electrons to the magnetic 
moment. We find 



A/z = 0.71 /J.B- 



(43) 



Since in our model the 4f moments have a fixed value of 
7 /iB the total moment amounts to 7.71 ^b very close 
to the experimental value of 7.63 /^g^- Our value is a 
bit smaller than that from the LDA+/7 calculation in 
Ref.El 

The procedure explained in the preceding sections al- 
lows for a determination of the full temperature depen- 
dence of the energy spectrum and the magnetic proper- 
ties of Gd. The selfconsistent evaluation yields a ferro- 
magnetic low-temperature phase with astonishing precise 
key-data. Fig. ^exhibits the magnetization curve in com- 
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FIG. 4: Magnetization curve for ferromagnetic Gadolinium as 
function of the temperature. The triangles represent experi- 
mental data taken from Ref. ,2. 



parison to the experimental data of Ref.|2- There is prac- 
tically an exact agreement. In particular the calculated 
Curie temperature, 



T c = 294.1 K, 



(44) 



being known as a very sensitive entity of magnetism, 
hardly deviates from the experimental value of 293.2 K. 
Note that there is in principle no fitting parameter in our 
theory, even the exchange constant J (I28|l is taken from 
the LDA input. We therefore have to conclude that the 

with the effective 



modified RKKY theory (Section llV B|) 

exchange integrals being functionals of the conduction 
electron selfenergy, describes the ferromagnetism of Gd 
in an absolutely convincing manner. 

Since we did not consider a direct exchange interac- 
tion between the localized 4f moments the induced spin 
polarization of the conduction electrons mediates the in- 
direct coupling. The "a priori" only slightly correlated 
5d/6s/6p band states therefore exhibit a distinct tem- 
perature dependence as can be seen for the total quasi- 
particle density of states in Fig. [SJ The T — splitting 
is responsible for the band contribution (43) to the to- 
tal magnetic moment. With increasing temperature the 
induced splitting reduces steadily collapsing at Tq. The 
shift is not at all rigid (" Stoner like" ) , but with clear de- 
formations. The latter point to a substantial influence 
of nonlinear effects such as magnon emission and absorp- 
tion and magnetic polaron formation, in particular what 
concerns the d states. The lower edge of the | spectrum, 
predominantly built up by 6s states (Fig. shows a 
red shift upon cooling below Tq as it is typical for local- 
moment systems, first observed for insulators and semi- 
conductors such as EuO and EuS^. The temperature 
behaviour at the chemical potential is not so clear. 
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FIG. 5: Quasiparticle density of states of the valence and 
conduction bands of Gd as function of the energy (E — /i) (/i: 
chemical potential) for four different temperatures. The total 
densities of states consist of 5d, 6s, and 6p contributions. 
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FIG. 6: Spin resolved single-electron spectral density of Gd 
at the r point as function of the energy (E — /i) for three 
different temperatures. Upper half: j spectrum, lower half: J, 
spectrum. 



The single-electron spectral density (jl 1|> represents the 
bare line shape of an angle- and spin-resolved photoc- 
mission experiment. Pronounced peaks in the spectral 
density define the quasiparticle band structure. For four 
high-symmetry points (F, A, H, M) we have calculated 
the energy dependence of the spectral density in the va- 
lence and conduction band region. The results for three 
different temperatures (T = 0, 200, 295 K) are repre- 
sented in Figs. El to El The T = 0-| spectra always con- 
sist of relatively sharp peaks pointing at quasiparticles 
with long, sometimes even infinite lifetimes. In case of 
infinite lifetime (real selfenergy) the spectral density is 
a (5-function. For plotting reasons we have then added 
a small imaginary part (iA; A = 0.01) to the electronic 
selfenergy. For empty energy bands the | spectrum would 
consist at T — exclusively of <5-peaks. This is just the 
exact limiting case discussed in Sec. MI Bl It means noth- 
ing else than that a | electron cannot undergo any scat- 
tering process if the localized 4f moments are ferromag- 
netically saturated. However, for finite and not fully spin- 
polarized band occupations some spin exchange processes 
may happen giving rise to slight quasiparticle dampings. 
Fig. shows the energy dependence of the spectral den- 
sity at the r point for three different temperatures. T = 
means ferromagnetic 4f saturation (Fig.0J while T = 295 
K is slightly above the calculated Curie temperature (|44H . 
At T = 200 K the moment system is partially ordered. 
The low energy peaks belong to 6s states (Figs.nand|21l. 
They are spin split in the ferromagnetic phase, where the 
induced exchange splitting diminishes continuously with 
increasing temperature, collapsing at T = Tq ("Stoner- 
like" behaviour). That agrees with the photoemission 
data of Kim et aln Similar temperature behaviour is 
found for the other quasiparticle peaks, too, and also 
for the other symmetry points A, H and M (Figs. [7] |8| 
EJ). These theoretical results contradict a bit our previ- 
ous investigation^ according to which in some cases a 
persisting splitting in the paramagnetic phase should be 



possible. The ambiguity comes along with the necessary 
decomposition of the total spectrum into non-degenerate 
subbands. That can be done, in principle, in different 
ways, and, at least in our opinion, it is not "a priori" 
clear which is the correct procedure. In this work we 
have used a method that retains the full atomic-orbital 
symmetry. The resulting rather broad subbands (Fig. EJ) 
cause correspondingly small effective exchange couplings 
J/W. The selfenergy M ka (E) (O is then dominated 
by the first term and therewith relatively close to the 
mean-field solution of the sf- (Kondo-lattice-) model. 



of 



A 

J 


iv 






l! J v 


T = 0K 

T = 200 K 

T - 295 K 






■ 







-6 -4 -2 2 4 6 8 10 



(E-H) [eV] 

FIG. 7: The same as in Fig. E]but for the A point. 

The " Stoner collapsing" is a typical feature of the weak- 
coupling (mean-field) region. The band decomposition 
used in Ref. |2f| leads to substantially smaller subband 
widths and therewith to stronger effective exchange cou- 
plings. 

A general observation is that with increasing tempera- 
ture the excitation peaks are getting broader, i.e. quasi- 
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FIG. 8: The same as in Fig. [S]but for the H point. 




FIG. 9: The same as in Fig.|^]but for the M point. 



particle lifetimes decrease. Raising temperature means 
enhancing magnon densities and therewith a higher prob- 
ability for electron-magnon spinflip scattering. The d-like 
states near and around the chemical potential /i exhibit 
stronger correlation effects than the low-energy 6s states 
or high-energy 6p states, again due to the larger effective 
coupling constant J/W. Contrary to the H and M points 
(Figs.|S10) there is no spectral weight at the chemical po- 
tential for the F and A points (Figs. El 0). 



many-body approach has previously been developed and 
tested in several model studies. It consists of a moment- 
conserving decoupling approach for the single-electron 
Green function, which fulfills a maximum number of ex- 
act limiting cases, and a modified RKKY theory for the 
localized moment system. The effective exchange inte- 
grals between the localized spins turn out to be func- 
tionals of the electronic selfenergy. In the weak coupling 
limit the approach agrees with the conventional RKKY 
theory. 

As single-electron (Bloch) energies we have used the 
results of an ASW band structure calculation therewith 
guaranteeing that all those interactions which are not 
explicitly covered by the Kondo-lattice model are taken 
into account in a rather realistic manner. An exact lim- 
iting case of the model could be exploited to avoid the 
well-known double counting problem. In a strict sense 
the method does not contain any really free parameter. 
The 4f-5d exchange coupling constant J, which enters 
the theory via the Kondo-lattice model, is fitted by the 
LDA input. 

The results of our theoretical investigation agree as- 
tonishingly well with the experimental data of Gd. The 
selfconsistent approach predicts correctly a ferromagnetic 
low-temperature phase. The magnetic T = moment 
is with 7.71 /iB very close to the experimental value of 
7.63 fiB- Even the extremely sensitive Curie temperature 
hardly deviates from the real Gd value (theory: 294.1 K, 
experiment: 293.2 K). The valence and conduction bands 
exhibit a remarkable induced temperature dependence. 
The T = exchange splitting explains the excess moment 
of 0.63 fiB (or 0.71 /zb), that cannot be ascribed to the 
seven 4f electrons. The temperature dependence of the 
exchange splitting roughly scales with the macroscopic 
magnetization collapsing at Tq ( "Stoner-behaviour" ) as 
it has been observed in photocmission experiments. Cor- 
relation effects lead to a distinct temperature dependence 
of the quasiparticle damping. 

We believe that the proposed combination of a careful 
many-body treatment of a proper theoretical model with 
an ab initio band structure calculation yields a rather re- 
alistic description of the ferromagnetic 4f metal Gadolin- 
ium. 
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